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Abstract 

This paper is speculated to propose a class of shrinkage estimators for shape parameter 
f in failure censored samples from two-parameter Weibull distribution when some ‘apriori’ or 
guessed interval containing the parameter f is available in addition to sample information and 
analyses their properties. Some estimators are generated from the proposed class and compared 
with the minimum mean squared error (MMSE) estimator. Numerical computations in terms of 
percent relative efficiency and absolute relative bias indicate that certain of these estimators 
substantially improve the MMSE estimator in some guessed interval of the parameter space of f, 
especially for censored samples with small sizes. Subsequently, a modified class of shrinkage 
estimators is proposed with its properties. 
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1. INTRODUCTION 

Identical rudiments subjected to identical environmental conditions will fail at different and 
unpredictable times. The ‘time of failure’ or ‘life length’ of a component, measured from some specified 
time until it fails, is represented by the continuous random variable X. One distribution that has been used 
extensively in recent years to deal with such problems of reliability and life-testing is the Weibull 
distribution introduced by Weibull(1939), who proposed it in connection with his studies on strength of 
material. 

The Weibull distribution includes the exponential and the Rayleigh distributions as special cases. 
The use of the distribution in reliability and quality control work was advocated by many authors following 
Weibull(1951), Lieblin and Zelen(1956), Kao(1958,1959), Berrettoni(1964) and Mann(1968 A). 
Weibull(1951) showed that the distribution is useful in describing the ‘wear-out’ or fatigue failures. 


Kao(1959) used it as a model for vacuum tube failures while Lieblin and Zelen(1956) used it as a model for 
ball bearing failures. Mann(1968 A) gives a variety of situations in which the distribution is used for other 
types of failure data. The distribution often becomes suitable where the conditions for “strict randomness” 
of the exponential distribution are not satisfied with the shape parameter # having a characteristic or 


predictable value depending upon the fundamental nature of the problem being considered. 


1.1 The Model 
Let x), x2, ..., X, be a random sample of size n from a two-parameter Weibull distribution, 


probability density function of which is given by : 
f(x;a,8) =pa vx exp{-(x/a)’}; x>0,a>0,8>0 
(1.1) 


where a@ being the characteristic life acts as a scale parameter and # is the shape parameter. 
The variable Y=1Inx follows an extreme value distribution, sometimes called the log-Weibull 


distribution [e.g. White(1969)], cumulative distribution function of which is given by : 


F(y)=1 ox oxo"), o<y<w0,-o<u<0,b>0 


(1.2) 
where b= 1/8 and u=In @ are respectively the scale and location parameters. 

The inferential procedures of the above model are quite complex. Mann(1967 A,B, 1968 B) 
suggested the generalised least squares estimator using the variances and covariances of the ordered 


observations for which tables are available up to n = 25 only. 


1.2 Classical Estimators 
Suppose x), X2, ..., Xm be the m smallest ordered observations in a sample of size n from Weibull 


distribution. Bain(1972) defined an unbiased estimator for 5 as 


h — y, BS 
ae i=l nk 
(m,n) 
(1.3) 
1 m-1 
where a= (4) #/ (0 rive ). 
(1.4) 


J; —u 


and Vv; = are ordered variables from the extreme value distribution with u=0O andb= 


AN 
1.The estimator 5, is found to have high relative efficiency for heavily censored cases. Contrary to this, 


AN 
the asymptotic relative efficiency of b,, is zero for complete samples. 


Engelhardt and Bain(1973) suggested a general form of the estimator as 


a ~ J; -y,| 
j,--5 | bovall 
j=1 | (g,m,n) 


(1.5) 


AN 
where g is a constant to be chosen so that the variance of bg is least and K(g mn) is an unbiasing constant. 
AN 


& 


The statistic has been shown to follow approximately y’ - distribution with h degrees of freedom, 


where h = 2) Var(b - /b) . Therefore, we have 
n— iP ®T(h/2)+j 
a” |-A( 75) eee 


= — j= 1,2 
B?\h-2 to 
(1.6) 
ae ey) ’ 2 i 
where {= a is an unbiased estimator of @ with Var ()= a ) and t = hb, having density 
1 B . ca (h/2)-1 
th=—T| = — |t ;t>0. 
ft) waa Hed 


The MMSE estimator of £, among the class of estimators of the form C @ ; C being a constant for 


which the mean square error (MSE) of Cf is minimum, is 


8 h-4 
M eo 
(1.7) 
having absolute relative bias and relative mean squared error as 
an. i 2) 
(1.8) 
and RMSE | By b=, 
h-2 


respectively. 


1.3 Shrinkage Technique of Estimation 

Considerable amount of work dealing with shrinkage estimation methods for the parameters of the 
Weibull distribution has been done since 1970. An experimenter involved in life-testing experiments 
becomes quite familiar with failure data and hence may often develop knowledge about some parameters of 
the distribution. In the case of Weibull distribution, for example, knowledge on the shape parameter / can 


be utilised to develop improved inference for the other parameters. Thompson(1968 A,B) considered the 


problem of shrinking an unbiased estimator € of the parameter € either towards a natural origin 5 or 
towards an_ interval (E, 6, and suggested the shrunken estimators A&+(1- h)s, and 


‘ +6§ 
hé nf) where 0 < h < 1 is a constant. The relevance of such type of shrunken 


estimators lies in the fact that, though perhaps they are biased, has smaller MSE than € for & in some 


gts 
interval around 5, or [Sc , as the case may be. This type of shrinkage estimation of the Weibull 


parameters has been discussed by various authors, including Singh and Bhatkulikar(1978), Pandey(1983), 
Pandey and Upadhyay(1985,1986) and Singh and Shukla(2000). For example, Singh and 
Bhatkulikar(1978) suggested performing a significance test of the validity of the prior value of @ (which 


they took as 1). Pandey(1983) also suggested a similar preliminary test shrunken estimator for £. 


In the present investigation, it is desired to estimate ( in the presence of a prior information 


available in the form of an interval (B, ‘5 B,) and the sample information contained in (2 . Consequently, 
this article is an attempt in the direction of obtaining an efficient class of shrunken estimators for the scale 
parameter 2. The properties of the suggested class of estimators are also discussed theoretically and 


empirically. The proposed class of shrunken estimators is furthermore modified with its properties. 


2. THE PROPOSED CLASS OF SHRINKAGE ESTIMATORS 


* 


Consider a class of estimators @ (eas for B in model (1.1) defined by 


Biase [Axe Ga ee | 


(2.1) 
where p and q are real numbers such that p#0 and gq > 0, w is a stochastic variable which may in 


* 


particular be a scalar, to be chosen such that MSE of £ tH) is minimum. 


Assuming w a scalar and using result (1.6), the MSE of Bi, 4) is given by 


wef = fad eae 2 ) Hee 


h—-2 T(h/2) 


+{qA- Hwa”( z J P(h/2) +22] 
h-2)  F(h/2) 


(2.2) 


where A= [EPs : 
2p 


Minimising (2.2) with respect to w and replacing £ by its unbiased estimator /, we get 


LP ae 


ae ep 
(Axe) 
2 
(2.3) 
_(h-2)” T[(n/2) +p] 
where wip)=[ 2 T[(A/2)+2p]’ 
(2.4) 


lies between 0 and 1, {i.e.,0 < w(p) < 1} provided gamma functions exist, i.c., p >(—A/2). 


Substituting (2.3) in (2.1) yields a class of shrinkage estimators for # in a more feasible form as 


Bae (+?) w(p) + os Ba (l-w(p)}. 


(2.5) 


2.1 Non-negativity 


Clearly, the proposed class of estimators (2.5) is the convex combination of {(h = 2)/ t } and 


{q(6, +B, )/2} and hence Box is always positive as {(h —2)/t} > 0 and q>0. 


2.2 Unbiasedness 
If w(p) = 1, the proposed class of shrinkage estimators B, p,q) turns into the unbiased estimator B, 


otherwise it is biased with 


Bis {B,,4)|=B {4d —1) wo) 26 


and thus the absolute relative bias is given by 


ARB {Beswf = | \gA on 1} [I — w(p)| |: 
(2.7) 


The condition for unbiasedness that w(p) = 1, holds iff, censored sample size m is indefinitely 
large, i.e., m — 00. Moreover, if the proposed class of estimators B (pq) tums into B then this case does not 
deal with the use of prior information. 


A more realistic condition for unbiasedness without damaging the basic structure of B pq) and 
utilises prior information intelligibly can be obtained by (2.7). The ARB of B toa is zero when g = Ae (or 


AG"): 


2.3 Relative Mean Squared Error 
The MSE of the suggested class of shrinkage estimators is derived as 


use [Pnof= P| Wawa? +2MPE) aw 


and relative mean square error is therefore given by 


° 2 ee : 
pe {Beruf ={gA-1'l—w(p)}° + aa , 


(2.9) 


It is obvious from (2.9) that RMSE Boro | is minimum when g = At (or A= qo ). 


2.4 Selection of the Scalar ‘p’ 
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The convex nature of the proposed statistic and the condition that gamma functions contained in 
w(p) exist, provides the criterion of choosing the scalar p. Therefore, the acceptable range of value of p is 
given by 
{p| 0<w(p)<1 and p>(-h/2)}, Vn.m. (2.10) 


2.5 Selection of the Scalar ‘q’ 

It is pointed out that at g = A’ , the proposed class of estimators is not only unbiased but renders 
maximum gain in efficiency, which is a remarkable property of the proposed class of estimators. Thus 
obtaining significant gain in efficiency as well as proportionately small magnitude of bias for fixed A or 
for fixed (B, , B) and (B > / B), one should choose g in the vicinity of g = A". Itis interesting to note 
that if one selects smaller values of g then higher values of A leads to a large gain in efficiency (along 
with appreciable smaller magnitude of bias) and vice-versa. This implies that for smaller values of q, the 
proposed class of estimators allows to choose the guessed interval much wider, i.e., even if the 


experimenter is less experienced the risk of estimation using the proposed class of estimators is not higher. 


This is legitimate for all values of p. 


2.3 Estimation of Average Departure: A Practical Way of selecting q 


The quantity A= {(B, +B,)/2B}, represents the average 


departure of natural origins f#, and #, from the true value 


f. But in practical situations it is hardly possible to get 
an idea about A. Consequently, an unbiased estimator of 
A is proposed, namely 


; ‘ (8, ‘#.)) T(h/2) 
A= ; 
4 T[(A/2) +1] 


(2:12) 


In section 2.5 it is investigated that, if Gah"; the 


suggested class of estimators yields favourable results. 
Keeping in view of this concept, one may select g as 
=k =| 4 {rie 
t(B,+B,)J T(h/2) 


Here this is fit for being quoted that this is the 


(2.13) 


criterion of selecting gq numerically and one should 
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carefully notice that this doesn’t mean gq is replaced by 


(213); ta fie 


3. COMPARISION OF ESTIMATORS AND EMPIRICAL STUDY 
James and Stein(1961) reported that minimum MSE is a highly desirable property and it is 
therefore used as a criterion to compare different estimators with each other. The condition under which the 


proposed class of estimators is more efficient than the MMSE estimator is given below. 
MSE 12. ms does not exceed the MSE of MMSE estimator / ,, if - 
(-VG)q" <A < (1+ VG)q" 
(3.1) 


2 1 {w(p)} 
{1-w(p)}° | (@-2)  (h-4) 


Besides minimum MSE criterion, minimum bias is also important and therefore should be 


where G= 


incorporated under study. Thus, ARB B pai) his less than ARB ‘B Fy, jit - 


2 2 
1 See rae | a 
| (h- 2)(I = Ga) ry — | i (A- 2y(I ~ Wop) i 


(3.2) 


3.1 The Best Range of Dominance of A 
The intersection of the ranges of A in (3.1) and (3.2) gives the best range of dominance of A 


denoted by A In this range, the proposed class of estimators is not only less biased than the MMSE 


Best * 


estimator but is more efficient than that. The four possible cases in this regard are: 


2 2 
(i) if i (h-2) i 7 | < ( JG) and t + (h-2) i = =| < ( + JG) then 


| 2 = 
A best {i 7 VG}4 ; i (h-2)[1- wall! | 


. 2 
(ii) if ; ool < ( 1G) and (+/G)< i 


is the same as defined in (3.1). 


+ 2 | then 
(h—2)[l- w(p)] 
A 


Best 
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2. 2 
ii. alee and (1+ /G ; | hen 
De | I< ial \\+ V6) < “GDI 


2 ay ioe i 
ne t al! A iGhe | 


2, 2 
iv if ae. 1 and 41 1 VG. hen 
wm at| )<| sal een aaron lid +¥6): 


is the same as defined in (3.2). 


A 


Best 


3.2 Percent Relative Efficiency 


To elucidate the performance of the proposed class of estimators on with the MMSE 


N N N 
estimator f ,, , the Percent Relative Efficiencies (PREs) of B,,,,. with respect to B,, have been computed 


by the formula: 
2(h —4) 
(h—2)[(gA-1)"{1—w(p)} h—4) + 2h)? 


PRE Baye By f= 


he 


A 


(3.5) | The PREs of Be 5 with respect to BY and ARBs of B., 4) for fixed n = 20 and different values 


of p,g,m A, (= B, /B) and A, (= B, /B) or A are compiled in Table 3.1 with corresponding values of / 


[which can be had from Engelhardt(1975)] and w(p). The first column in every m corresponds to PREs and 


the second one corresponds to ARBs of une The last two rows of each set of g includes the range of 


dominance of A and A,,., 


. The ARBs of BY has also been given at the end of each set of table. 
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Table 3.1 


PREs of proposed estimator B oa with respect to MMSE estimator B,, and ARBs of B coef 


Range of A> 


Agest > 


0.50 
1.00 
1.50 
2.00 
2.50 
3.00 
3.50 
4.00 


0.15 
0.50 
1.00 
1.50 
2.00 
2.50 
3.00 
3.50 
4.00 
Range of A> 


ABest > (0.97, 1.70) (1.01, 1.66) (1.03, 1.64) (1.04, 1.63) 
ARB of MMSE Estimator > 0.2259 0.1463 0.1061 0.0820 
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Table 3.1 continued ... 


3.50 
4.00 


0.15 
0.50 
1.00 
1.50 
2.00 
2.50 
3.00 
3.50 
4.00 


0.15 
0.50 
1.00 
1.50 
2.00 
2.50 
3.00 
3.50 
4.00 
Range of A> 


Apest > (0.00, 2. 67) 0. 00, 2. 67) 0. 00, 2. 67) (0.00, 2. 67) 


ARB of MMSE Estimator > 0.2259 0.1463 0.1061 0.0820 
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Table 3.1 continued ... 


0.50 
1.00 
1.50 
2.00 
2.50 
3.00 
3.50 
4.00 
Range of A> 


Ais (0.37, 2.30) (0.47, 2.20) (0.52, 2.14) (0.44, 2.23) 


ARB of MMSE Estimator > 0.2259 0.1463 0.1061 0.0820 
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Table 3.1 continued ... 


0.50 
1.00 
1.50 
2.00 
2.50 
3.00 
3.50 
4.00 
Range of A> 


Apest > (0.89, 1.77) (0.99, 1.68) (1.03, 1.64) (0.99, 1.68) 


ARB of MMSE Estimator > 0.2259 0.1463 0.1061 0.0820 
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3000.00 


2500.00 n= 20 
m = 6 
~ Sea 
S 2000.00 q = 0.25 
a 
S = 
= 1500.00 | | —#— PRE 
: a 7 —e— ARB*1000 
a 1000.00 / \ seeeee ARB(MMSE Esti.) 
/ : PRE Cut-off Point 


0.05 1 2 3 4 5 6 7 8 


It has been observed from Table 3.1, that on keeping m, p, q fixed, the relative efficiencies of the 
proposed class of shrinkage estimators increases up to A = q'', attains its maximum at this point and then 
decreases symmetrically in magnitude, as A increases in its range of dominance for all n, p and q. On the 
other hand, the ARBs of the proposed class of estimators decreases up to A = q'', the estimator becomes 
unbiased at this point and then ARBs increases symmetrically in magnitude, as A increases in its range of 
dominance. Thus it is interesting to note that, at ¢ = A”', the proposed class of estimators is unbiased with 
largest efficiency and hence in the vicinity of g = A"! also, the proposed class not only renders the massive 
gain in efficiency but also it is marginally biased in comparison of MMSE estimator. This implies that q 


plays an important role in the proposed class of estimators. The following figure illustrates the discussion. 


Figure 3.1 


The effect of change in censored sample size m is also a matter of great interest. For fixed p, g and 
A, the gain in relative efficiency diminishes, and ARB also decreases, with increment in m. Moreover, it 
appears that to get better estimators in the class, the value of w(p) should be as small as possible in the 
interval (0,1]. Thus, to choose p one should not consider the smaller values of w(p) in isolation, but also the 


wider length of the interval of A. 
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4. MODIFIED CLASS OF SHRINKAGE ESTIMATORS AND ITS PROPERTIES 
The proposed class of estimators B, pq) 18 not uniformly better than B . It will be better if B, and 
B., are in the vicinity of true value fi . Thus, the centre of the guessed interval (B, +B, \/ 2 is of much 


importance in this case. If we partially violate this, i.e., only the centre of the guessed interval is not of 


much importance, but the end points of the interval B, and {, are itself equally important then we can 


propose a new class of shrinkage estimators for the shape parameter [i by using the suggested class 


Bas ee 
B, , if ¢>[(A-2)/B,] 

Bann =a) tv) +f PPE ton}. if [in-29/Ba]s 4 < [-29/)] 
B, , if t<[(h-2)/B,] 

(4.1) 

which has 


Biff =A aft A(8 be ord 7 i{na5 i} 
+attwtoi(mk}- ang bas (noe) 


(4.2) 
and 
MSEf,,.. |= ae —1)’ -A,(A, 2) i[niek} +0208; 2) i{na.2) 
o(h-2 h h 
+{w(p)} Gite 2| i[n.5 2} 
+adti—woo n.)-1[ nae] Heat o(o)}-2 
h h 
raeo)4t( nd 7 i[ns4 i} Haatt—wen}—) 
(4.3) 
where -(3-1)a: n -(3 a3 and I(n,@)= : feud 
1 1? 2 5 2 ] T(@) , : 


19 


This modified class of shrinkage estimators is proposed in accordance with Rao(1973) and it 
seems to be more realistic than the previous one as it deals with the case where the whole interval is taken 


as apriori information. 


5. NUMERICAL ILLUSTRATIONS 
The percent relative efficiency of the proposed estimator B, pq) With respect to MMSE 


N 
estimator B,, has been defined as 


ee Bia By, = se a} x 100 
(p.9) 


(5.1) 
and it is obtained for n = 20 and different values of p, g, m, A, and A, (or A). The findings are 


summarised in Table 5.1 with corresponding values of h and w(p). 


Table 5.1 


PREs of proposed estimator B, aa with respect to MMSE estimator B 


n=20 


6 12 6 8 10 12 
10.8519]15.6740 
0.8537 0.7737 


50.80) 41.39 40.10 
117.60} 81.01 79.57 
261.72) 227.42 191.97 
548.60} 426.98 355.31 
649.95} 470.44 504.49 
268.31) 189.82 210.91 
80.46} 53.66 55.10 
50.84 
120.81 
298.17 
642.86 
626.09 
247.90 
78.41 

0.2 | 0.3 0.25 50.89) 41.24) 34.60} 30.19] 49.97} 39.95) 34.23) 30.59 

0.4 | 0.6 0.50 124.02} 83.01) 68.50} 63.81) 122.74) 82.68) 67.32} 62.50 

0.6 | 0.9 0.75 339.92} 282.24) 242.46) 197.73} 325.66) 266.36) 229.58] 192.68 
0.75 | 0.8 | 1.2 1.00 723.50} 510.42) 389.34) 316.87] 710.96) 504.67) 388.35} 317.53 

1.0 ) 1.5 125 566.19} 392.47) 312.16) 263.77] 597.64) 421.61] 337.17) 278.26 


20 


1.2 | 1.8 1.50 224.67) 161.95) 131.14} 111.81) 233.41) 169.19} 136.65} 114.63 
1.5 | 2.0 1.75 76.05) 51.59} 38.85) 30.83] 76.93} 52.14) 39.17) 30.95 
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Table 5.1 continued ... 


6 12 6 8 10 12 
10.8519]15.6740 
0.8537 0.7737 


34.18 34.81 
72.59 73.36 
95.44 102.24 
131.16 145.23 
243.10 273.81 
211.32 220.57 
57.89 57.95 


It has been observed from Table 5.1 that likewise Bias the PRE of Bess, with respect to B 


m 
decreases as censoring fraction (m/n) increases. For fixed m, p and q the relative efficiency increases up to 
a certain point of A, procures its maximum at this point and then starts decreasing as A increases. It 
seems from the expression in (4.3) that the point of maximum efficiency may be a point where either any 


one of the following holds or any two of the following holds or all the following three holds- 


(i) the lower end point of the guessed interval, i.e., 8, coincides exactly with the true value B, ice., 
A,=1. 

(ii) the upper end point of the guessed interval, i.e.,., departs exactly two times from the true value 
Bie, A,=2. 

ji) Az=q" 


This leads to say that on contrary to B (pq) there is much importance of A, and A, in addition to A. 


The discussion is also supported by the illustrations in Table 5.1. As well, the range of dominance of 
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average departure A is smaller than that is obtained for B tans but this does not humiliate the merit of 


B (pg) because still the range of dominance of A is enough wider. 


6. CONCLUSION AND RECOMMENDATIONS 

It has been seen that the suggested classes of shrunken estimators have considerable gain in 
efficiency for a number of choices of scalars comprehend in it, particularly for heavily censored samples, 
i.e., for small m. Even for buoyantly censored samples, i.e., for large m, so far as the proper selection of 
scalars is concerned, some of the estimators from the suggested classes of shrinkage estimators are more 
efficient than the MMSE estimators subject to certain conditions. Accordingly, even if the experimenter has 


less confidence in the guessed interval (B, B sy of B, the efficiency of the suggested classes of shrinkage 


estimators can be increased considerably by choosing the scalars p and g appropriately. 
While dealing with the suggested class of shrunken estimators By pq) it is recommended that one 


should not consider the substantial gain in efficiency in isolation, but also the wider range of dominance of 
A ,, because enough flexible range of dominance of A will leads to increase the possibility of getting 
better estimators from the proposed class. Thus it is recommended to use the proposed class of shrunken 


estimators in practice. 


REFERENCES 


BAIN, L. J. (1972) : Inferences based on Censored Sampling from the Weibull or Extreme-value 
distribution, Technometrics, 14, 693-703. 


BERRETTONI, J. N. (1964) : Practical Applications of the Weibull distribution, /ndustrial Quality 
Control, 21, 71-79. 


ENGELHARDT, M. and BAIN, L. J. (1973) : Some Complete and Censored Sampling Results for the 


Weibull or Extreme-value distribution, Technometrics, 15, 541-549. 


ENGELHARDT, M. (1975) : On Simple Estimation of the Parameters of the Weibull or Extreme-value 
distribution, Technometrics, 17, 369-374. 


JAMES, W. and STEIN, C. (1961) : (A basic paper on Stein-type estimators), Proceedings of the 4" 
Berkeley Symposium on Mathematical Statistics, Vol. 1, University of California Press, Berkeley, CA, 
361-379. 


KAO, J. H. K. (1958) : Computer Methods for estimating Weibull parameters in Reliability Studies, 
Transactions of IRE-Reliability and Quality Control, 13, 15-22. 


23 


KAO, J. H. K. (1959) : A Graphical Estimation of Mixed Weibull parameters in Life-testing Electron 
Tubes, Technometrics, 1, 389-407. 


LIEBLEIN, J. and ZELEN, M. (1956) : Statistical Investigation of the Fatigue Life of Deep Groove Ball 
Bearings, Journal of Res. Natl. Bur. Std., 57, 273-315. 


MANN, N. R. (1967 A) : Results on Location and Scale Parameter Estimation with Application to the 
Extreme-value distribution, Aerospace Research Labs, Wright Patterson AFB, AD.653575, ARL-67-0023. 


MANN, N. R. (1967 B) : Table for obtaining Best Linear Invariant estimates of parameters of Weibull 
distribution, Technometrics, 9, 629-645. 


MANN, N. R. (1968 A) : Results on Statistical Estimation and Hypothesis Testing with Application to the 
Weibull and Extreme Value Distribution, Aerospace Research Laboratories, Wright-Patterson Air Force 


Base, Ohio. 


MANN, N. R. (1968 B) : Point and Interval Estimation for the Two-parameter Weibull and Extreme-value 
distribution, Technometrics, 10, 231-256. 


PANDEY, M. (1983) : Shrunken estimators of Weibull shape parameters in censored samples, JEEE Trans. 
Reliability, R-32, 200-203. 


PANDEY, M. and UPADHYAY, S. K. (1985) : Bayesian Shrinkage estimation of Weibull parameters, 
IEEE Transactions on Reliability, R-34, 491-494. 


PANDEY, M. and UPADHYAY, S. K. (1986) : Selection based on modified Likelihood Ratio and 


Adaptive estimation from a Censored Sample, Jour. Indian Statist. Association, 24, 43-52. 
RAO, C. R. (1973) : Linear Statistical Inference and its Applications, Sankhya, Ser. B, 39, 382-393. 


SINGH, H. P. and SHUKLA, S. K. (2000) : Estimation in the Two-parameter Weibull distribution with 
Prior Information, JAPOR Transactions, 25, 2, 107-118. 


SINGH, J. and BHATKULIKAR, S. G. (1978) :Shrunken estimation in Weibull distribution, Sankhya, 39, 
382-393. 


THOMPSON, J. R. (1968 A) : Some Shrinkage Techniques for Estimating the Mean, The Journal of 


American Statistical Association, 63, 113-123. 


THOMPSON, J. R. (1968 B) : Accuracy borrowing in the Estimation of the Mean by Shrinkage to an 


Interval , The Journal of American Statistical Association, 63, 953-963. 


WEIBULL, W. (1939) : The phenomenon of Rupture in Solids, Ingenior Vetenskaps Akademiens 
Handlingar, 153,2. 


WEIBULL, W. (1951) : A Statistical distribution function of wide Applicability, Journal of Applied 
Mechanics, 18, 293-297. 


24 


WHITE, J. S. (1969) : The moments of log-Weibull order Statistics, Technometrics,11, 373-386. 


25 


